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ARTICLE INFO ABSTRACT 


Keywords: Porcine deltacoronavirus (PDCoV) is a novel porcine enteric coronavirus that causes diarrhea, vomiting and 
Porcine deltacoronavirus dehydration in piglets. This newly virus has spread rapidly and has caused serious economic losses for pig 
Prevalence 


industry since the outbreak in USA in 2014. In this study, 430 faecal and intestinal samples (143 faecal samples 
and 287 intestinal samples) were collected from individual pigs with diarrhea and 211 serum samples were also 
collected from the sows with mild diarrhea in 17 regions in Henan province, China from April 2015 to March 
2018. The RT-PCR detection indicated that the infection of PDCoV was high up to 23.49% (101/430), and co- 
infection with PEDV were common (60.40%, 61/101) in Henan pigs. The prevalence of PDCoV in suckling 
piglets was the highest (36.43%, 94/258). We also found that PDCoV could be detected in sows faeces and sera 
while the sows showed mild, self-limited diarrhea in clinic. The complete genomes of 4 PDCoV Henan strains 
(CH-01, HNZK-02, HNZK-04, HNZK-06) were sequenced and analyzed. Phylogenetic analysis based on the 
complete genome, spike and nucleocapsid gene sequences revealed that the PDCoV Henan strains were closely 
related to other PDCoV reference strains that located in the Chinese clade. Furthermore, the phylogenetic 
analysis showed PDCoV CH-01 strain was closely related to CHN-HB-2014 strain and HKU15-44 strain, while the 
other PDCoV Henan strains were more related to PDCoV CHJXNI2 and CH-SXD1-2015 strains, indicating that 
the ancestor of these sequenced strains may different. These results would support the understanding of the 
prevalence and evolution characteristics of PDCoV in China. 


Complete genome 
Phylogenetic analysis 
Evolution 


1. Introduction 


Coronaviruses are single-stranded positive-sense, enveloped RNA 
viruses that can be classified into four genera, Alphacoronavirus, 
Betacoronavirus, Gammacoronavirus and Deltacoronavirus (Ma et al., 
2015). Coronaviruses from different genera can infect diverse host 
species, mainly including mammals and birds (Woo et al., 2010; Chan 
et al., 2013). Deltacoronavirus, a novel genus, was first found in Asian 
leopard cat and Chinese ferret badgers in Hong Kong and then identi- 
fied in a variety of avian species (Dong et al., 2007). In 2012, the 
porcine deltacoronavirus (PDCoV) was first discovered in pigs in Hong 
Kong (Woo et al., 2012). To date, six coronaviruses have been identified 
from pigs, including the Alphacoronavirus porcine epidemic diarrhea 
virus (PEDV), transmissible gastroenteritis virus (TGEV), porcine re- 
spiratory coronavirus (PRCV) and swine acute diarrhea syndrome cor- 
onavirus (SADS-CoV), the  Betacoronavirus genus porcine 


haemagglutinating encephalomyelitis coronavirus (PHEV), and the 
Deltacoronavirus genus PDCoV (Jung et al., 2016b; Gong et al., 2017; 
Pan et al., 2017; Hu et al., 2018; Zhou et al., 2018). The clinical 
symptom caused by PDCoV in piglets is similar to that of the PEDV and 
TGEV infection, including severe diarrhea, vomiting and dehydration 
(Hu et al., 2016; Jung et al., 2016a). However, the clinical symptom 
and mortality caused by PDCoV are milder than PEDV infection in 
piglets (Jung et al., 2016a). Furthermore, PDCoV commonly co-infected 
with PEDV and TGEV in natural infections and has seriously harmed the 
development of the pig industry (Hu et al., 2015). Currently, there are 
no effective reagents and vaccines to control PDCoV infection. 

PDCoV has the smallest genome size among coronaviruses and the 
full-length genome is approximately 25.4kb with a 5’ cap and a 3’ 
polyadenylated tail (Fang et al., 2016). The two-thirds of PDCoV 
genome at the 5’ end encodes genome replication proteins (ORFlab) 
and one-third at the 3’ end encodes four mainly structural proteins, 
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spike glycoprotein (S), envelope protein (E), membrane protein (M), 
and nucleocapsid protein (N), and three accessory proteins, including 
nonstructural protein 6 (NS6), NS7, and NS7a (Woo et al., 2012; Jang 
et al., 2017). The PDCoV S protein is a glycoprotein on the viral surface 
and plays an essential function in mediating viral entry into host cells 
by binding with host receptors (Graham and Baric, 2010; Shang et al., 
2017). In addition, PDCoV S gene encodes the most antigenic protein 
and frequently used for studying the genetic relationships among dif- 
ferent PDCoV strains and the epidemiological investigation (Shang 
et al., 2017). The highly conserved coronaviruses N protein produces 
largely in infected cells and plays a vital function in viral replication 
and pathogenesis. The recent study showed that the multifunctional 
role of PDCoV N protein is similar to other coronaviruses (Lee and Lee, 
2015). 

In February 2014, PDCoV was first detected in swine herds with 
diarrhea in Ohio, USA (Wang et al., 2014a). This novel coronavirus was 
then further detected in more than 20 US states (Marthaler et al., 2014; 
Wang et al., 2014b; Sinha et al., 2015). Subsequently, PDCoVs were 
detected in pigs in Canada, South Korea, mainland China, Thailand and 
Laos (Lee and Lee, 2014; Janetanakit et al., 2016; Saeng-Chuto et al., 
2017; Ajayi et al., 2018). In China, PDCoV infection in pig farms with 
diarrhea occurred in some provinces and regions since 2014 (Wang 
et al., 2015; Dong et al., 2016; Zhai et al., 2016). Previous study also 
showed that PDCoV was detected in diarrheic pig samples collected as 
early as 2004 in China (Dong et al., 2015). Sequence analyses suggested 
that all the known PDCoVs strains have a high nucleotide similarity 
with the earliest discovered strains, HKU15-44 and HKU15-155 (Jung 
et al., 2016b), and all reported China PDCoV strains are closely related 
with each other and located in the same clade (Song et al., 2015; Liu 
et al., 2018; Mai et al., 2018; Wang et al., 2018). Henan province is a 
major pig-raising region in China. According to the data of the National 
Bureau of Statistics of China, the number of pigs in stock in the first half 
year of 2018 in Henan province was about 42 million with the per- 
centage of 10.2% in the whole stock pigs in China. However, little is 
reported regarding the prevalence and epidemiology of PDCoV in 
Henan province. 

In this study, we aimed to investigate the PDCoV prevalence in the 
diarrheal pigs and full-length genome sequence properties of PDCoV in 
Henan province. 430 faecal and intestinal samples were collected from 
pigs with diarrhea and vomiting in 17 regions for PDCoV testing. 211 
serum samples were also collected from the sows with mild diarrhea to 
detect PDCoV infection. The characteristic of PDCoV prevalence in the 
diarrheal pigs in different regions of Henan province was analyzed. 
Then the complete genomes of 4 representative PDCoV Henan strains 
were sequenced, and the phylogenetic and evolutionary characteriza- 
tion of these four PDCoV strains were further analyzed. Our study will 
provide more information about prevalence and genetic evolution of 
PDCoV. 


2. Materials and methods 
2.1. Clinical sample collection 


From April 2015 to March 2018, 430 faecal and intestinal samples 
(143 faecal samples and 287 intestinal samples) were collected from 
pigs with diarrhea outbreaks that diagnosed by pig farm veterinarian in 
17 regions in Henan province, China. The number of sampled farms and 
pigs in different years and regions was as shown in Table 1. 211 serum 
samples were also collected from the sows with mild diarrhea in this 
period. The swine faeces and intestinal contents were diluted 10-fold 
with phosphate-buffered saline (PBS). The suspensions were vortexed 
and centrifuged at 1847 g at 4°C for 10 min. The clarified supernatants 
were collected and stored at — 80°C for RNA extraction. The research 
protocol for animal experiments of live pigs in this study was approved 
by the Animal Care and Use Committee of Henan Agricultural Uni- 
versity (Zhengzhou, China) and was performed in accordance with the 
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“Guidelines for Experimental Animals” of the Ministry of Science and 
Technology (Beijing, China). 


2.2. RNA extraction and PDCoV detection 


Total RNAs were extracted from samples using Trizol reagent 
(Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s in- 
structions. Synthesis of cDNA from the extracted RNA was carried out 
by using the Vazyme Reverse Transcription Kit (Nanjing, China). To 
detect PDCoV from the samples’ RNA, a pair of primers targeting the S 
gene of PDCoV (5’- GACCCTAAATCTGCCGTTAGAG-3’ and 5’-TGTTG 
GAGAGGTGAATGCTATG-3’) were designed based on the sequence of 
HKU15-155 strain (GenBank accession no. JQ065043), and the ex- 
pected size of the amplicon is 543 bp. PCR was performed in a 25 uL 
volume containing 1 uL cDNA, 12.5uL 2x Rapid Taq Master Mix, and 
0.5 pL of each primer (25 pmol/L). The reactions were subjected to 95°C 
for 5 min, followed by 35 cycles of 95°C for 50s, 58°C for 45 s, and 72°C 
for 30s, with a final extension step at 72°C for 10 min. The PCR pro- 
ducts were purified using the SanPrep Column DNA Gel Extraction Kit 
(Sangon Biotech, Shanghai, China). Both strands of the PCR products 
were sequenced twice with an ABI Prism 3700 DNA analyzer (Sangon 
Biotech, Shanghai, China) by using the forward and reverse PCR pri- 
mers. The sequences of the PCR products were compared with the 
known S gene sequences of PDCoVs in the GenBank database. The 
PDCoV positive test samples were then further used to detect PEDV and 
TGEV. To detect PEDV and TGEV from the samples’ RNA, two pair of 
primers targeting the M gene of PEDV (5’-TGGTGGTCTTTCAATC 
CTG-3’ and 5’-CTGAGTAGTCGCCGTGTTT-3’) and the N gene of TGEV ( 
5’-CGCTATCGCATGGTGAAG-3’ and 5’-GGATTGTTGCCTGCCTCT-3’) 
were designed based on the sequence of the JXNC1302 strain (GenBank 
accession no. KJ526109) and the H155 strain (GenBank accession no. 
GQ374566) respectively. The expected size of the amplicons are 298 bp 
for the M gene of PEDV and 324 bp for the N gene of TGEV. The PCR 
amplifications and products analysis were the same as that of PDCoV. 


2.3. Complete genome sequence analysis 


Four representative PDCoV Henan strains were selected from four 
different pig farms that the pigs had severe diarrhea at different times to 
sequence the complete genome by the Sanger method. These re- 
presentative PDCoV strains were named as CH-01, HNZK-02, HNZK-04 
and HNZK-06. To determine the full-length genomic sequence of the 
isolated PDCoV strains, 13 pairs of primers were designed based on the 
HKU15-155 strain (GenBank accession no. JQ065043) (Table 2). Viral 
RNA was extracted according to the aforementioned method. Thirteen 
overlapping fragments covering the entire viral genome were amplified 
on the conditions of 95°C for 5 min, followed by 35 cycles of 95 °C for 
1 min, 55°C for 45s, 72°C for 2 min and finally 72°C for 10 min. The 
PCR products were purified and cloned into pMD18-T vectors (TaKaRa, 
Dalian, China) following the manufacturer’s protocol. The positive 
clone of each amplicon were sequenced twice (Sangon Biotech, 
Shanghai, China). 


2.4. Sequence alignment and phylogenetic analyses 


Complete sequences of identified PDCoV strains in this study were 
obtained by assembling overlapping contigs followed by trimming off 
primer sequence, and then used in sequence alignment and phyloge- 
netic analyses with other PDCoV reference strains (the representative 
strains of USA, China, South Korea, Thailand, Lao and Vietnam). The 
information of the four isolate strains and other reference strains of 
PDCoV are shown in Table 3. The nucleotide and deduced amino acid 
sequences were aligned and analyzed using the Clustal W program of 
DNAStar 7.0 green (DNAstar, Madison, WI). Phylogenetic trees were 
constructed using the neighbor-joining method in MEGA 6.0 software 
(http://www.megasoftware.net/) with bootstrap analysis of 1000 


H. Zhang, et al. 


Table 1 
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The number of sampled farms and pigs in different years and regions in Henan province, China. 


Regions The number of sampled farms The number of sampled pigs 
2015 2016 2017 2018 Total 2015 2016 2017 2018 Total 
An yang 0 1 1 0 2 0 2 4 0 6 
He Bi 1 4 3 1 9 3 8 12 3 26 
Ji Yuan 0 0 2 0 2 0 0 6 0 6 
Kai Feng 4 2 5 1 12 15 6 14 2 37 
Luo Yang 1 5 1 0 7 1 17 6 0 24 
Luo He 0 0 1 0 1 0 3 0 3 
Nan Yang 3 3 4 1 11 4 20 17 2 43 
Ping Dingshan 3 1 3 1 8 7 6 18 6 37 
Pu Yang 2 2 6 1 11 10 16 51 2 79 
San Menxia 0 1 3 0 4 0 1 9 0 10 
Zheng Zhou 1 0 4 1 6 2 0 10 3 15 
Xin Xiang 1 3 1 0 5 2 7 5 0 14 
Xin yang 0 0 3 2 5 0 0 3 3 6 
Zhou Kou 2 3 4 2 11 5 8 27 31 71 
Zhu Madian 1 0 5 0 6 7 0 21 0 28 
Shang Qiu 1 1 4 1 vb 3 1 17 1 22 
Xu Chang 0 0 1 1 2 0 0 1 2 3 
Total 20 26 51 12 109 59 92 224 55 430 
replicates. PDCoV detection. The results showed that 101 PDCoV positive samples 


2.5. Nucleotide sequence accession numbers 


The complete genome sequences of PDCoV CH-01, HNZK-02, HNZK- 
04 and HNZK-06 strains were deposited in GenBank under accession 
numbers KX443143, MH708123, MH708124 and MH708125, respec- 


tively. 


3. Results 


3.1. Prevalence of PDCoV in clinical samples from diarrhoeic pigs 


In total, 430 porcine clinical samples with diarrhea were collected 
from 17 regions in Henan province from April 2015 to March 2018 for 


Table 2 


Primers used for amplification of the complete genome of PDCoV. 


were identified from the 430 suspected samples and the total positive 
proportion was 23.49% (101/430) as shown in Table 4. In addition, 
PDCoV could be detected in sows’ sera samples with the positive pro- 
portion of 25.59%. The positive proportion of PDCoV in different re- 
gions since April 2015 was shown in Fig. 1. PDCoV was detected out in 
14 regions out of the 17 regions, including the regions of Xuchang 
(100%), Luohe (100%), Zhengzhou (66.67%) and Zhoukou (42.25%), 
implying this novel coronavirus existed widely in pig farms in Henan 
province (Fig. 1). 

The positive proportions of PDCoV in different years and seasons 
were analyzed and shown in Fig. 2A-B. The results indicated that the 
prevalence of PDCoV in the diarrheal pigs was increasing from 2015 
(1.69%) to 2018 (81.82%). The PDCoV positive proportion in spring 
(20.57%) and winter (54%) was obviously higher than that in summer 


Primers Sequence (5’—3’) Nucleotide position” Fragment size (bp) 

PDC-1-F ACATGGGGACTAAAGATAAAAATTATAGC 1-29 3300 
PDC-3300-R GCTCATCGCCTACATCAGTG 3281-3300 

PDC-3091-F CGGATTTAAAACCACAGACT 3091-3110 1770 
PDC-4860-R ACGACTTTACGAGGATGAAT 4841-4860 

PDC-4741-F CTCCTGTACAGGCCTTACAA 4741-4760 1680 
PDC-6420-R TCACACGTATAGCCTGCTGA 6401-6420 

PDC-6291-F CTCAATGCAGAAGACCAGTC 6291-6310 1751 
PDC-8041-R CAGCTTGGTCTTAAGACTCT 8041-8060 

PDC-7920-F GGTACTGCTTCTGATAAGGAT 7920-7940 1741 
PDC-9660-R TAGGTACAGTTGTGAACCGA 9641-9660 

PDC-9463-F TACTCTTCACAGCCTTCAC 9463-9482 1117 
PDC-10579-R TTCTTCACAATAGTGTAGCG 105601-10579 

PDC-10560-F CGCTACACTATTGTGAAGAA 10560-10579 2288 
PDC-12847-R TTCGTAGGGCTCAAGATA 12829-12847 

PDC-12301-F TCCAGATGACCGTTGCGTAT 12301-12320 2652 
PDC-14952-R CCAACAGAGTCGGGTAAA 14935-14952 

PDC-14703-F CACCCATAACGAAGAACC 14703-14720 2483 
PDC-17185-R CCGATGAGTGTCGTAGCG 17168-17185 

PDC-17165-F TGCCGCTACGACACTCAT 17165-17182 2039 
PDC-19203-R TCCGCTAAAGGAGAATAGGTTG 19182-19203 

PDC-18485-F TGCTACCCAATCTTACAGT 18485-18503 2594 
PDC-21079-R GCAAATACTCCGTCTGAAC 21061-21079 

PDC-20761-F GTCTTACCGTGTGAAACCCC 20761-20780 3064 
PDC-23825-R AACCACGAGACTGTAAGCAA 23806-23825 

PDC-23804-F TTTTGCTTACAGTCTCGTGGTT 23804-23825 1629 
Oligo(dT) TITTTTTTTTTTTTTTTIT 25432 


* Nucleotide position of each primer was numbered based on the HKU15-155 strain (GenBank accession no. JQ065043). 
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Table 3 
The isolate strains and other reference strains of PDCoV used for sequence 
alignment and phylogenetic analysis. 


Strain Accession no. Country Isolated time 
CH-01 KX443143 China/Hebi 2016 
HNZK-02 MH708123 China/Zhoukou 2018 
HNZK-04 MH708124 China/Zhoukou 2018 
HNZK-06 MH708125 China/Zhoukou 2018 
CHN-AH-2004 KP757890 China/Anhui 2004 
CHJXNI2 KR131621 China/Jiangxi 2015 
CHN-HB-2014 KP757891 China/Hebei 2014 
CHN-JS-2014 KP757892 China/Jiangsu 2014 
CHN-HN-2014 KT336560 China/Henan 2014 
CH/JXJGSO1/2016 KY293677 China/Jiangxi 2016 
CH/SXD1/2015 KT021234 China 2015 
HKU15-155 JQ065043 China/Hong 2010 
Kong 
HKU15-44 JQ065042 China/Hong 2009 
Kong 
Ohio137 KJ601780 USA 2014 
Illinois134 KJ601778 USA 2014 
Illinois121 KJ481931 USA 2014 
Indiana453 KR265851 USA 2014 
Minnesotal59 KR265859 USA 2014 
Minnesota292 KR265864 USA 2014 
Minnesota KR265853 USA 2013 
NorthCarolina452 KR265858 USA 2014 
Arkansas61 KR150443 USA 2015 
Michigan/8977 KM012168 USA 2014 
Nebraska209 KR265860 USA 2014 
OhioCVM1 KJ769231 USA 2014 
NE3579 KJ584359 USA 2014 
USA/IA/2014 KJ567050 USA 2014 
PA3148 KJ584358 USA 2014 
MI6148 KJ620016 USA 2014 
OH1987 KJ462462 USA 2014 
SD3424 KJ584356 USA 2014 
Swine/Thailand/S5025 KU051656 Thailand 2015 
PDCoV/2016/Lao KX118627 Lao 2016 
Swine/Vietnam/Binh21/ KX834352 Vietnam 2015 
2015 
KNU16-07 KY364365.1 South Korea 2014 
KNU14-04 KM820765 South Korea 2014 
Table 4 


The detection result of PDCoV epidemic in different years and seasons in Henan 
province, China. 


Seasons Sample number Positive number of PDCoV 
2015 59 1 
Spring 27 0 
Summer 13 0 
Autumn 14 0 
Winter 16 3 
2016 92 22 
Spring 34 11 
Summer 18 4 
Autumn 15 0 
Winter 29 14 
2017 224 33 
Spring 132 19 
Summer 32 0 
Autumn 29 0 
Winter 55 37 
2018 55 45 
Spring 16 13 
Total 430 101 


(6.35%) and autumn (0%) (Fig. 2B). Moreover, the positive proportion 
of PDCoV infection in suckling piglets was 36.43%, which was the 
highest at different stages of pigsin Henan province. PDCoV infection 
was relatively low in nursery pigs with the positive proportion of 
7.14%. PDCoV positive proportion was the lowest in sows faeces 
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(2.59%) in our current study (Fig. 2C). Our detection also found that co- 
infections of PDCoV with PEDV were common in Henan diarrhoeic pigs. 
61 out of 101 PDCoV positive samples were proved to exist PEDV, with 
the positive proportion of 60.40%. In addition, TGEV was identified in 
6 out of 101 (5.94%) PDCovV positive samples (Fig. 2D). These data also 
suggested that the co-infections of PDCoV with PEDV and TGEV often 
occurred and had caused great damage to local pig farms. 


3.2. Complete genomic characterization of PDCoV strains in Henan, China 


The complete genome sequence of four representative PDCoV 
strains (CH-01, HNZK-02, HNZK-04, HNZK-06) were determined using 
the designed primers. The results showed that the entire genome of 
PDCoV CH-01 strain comprised 25,404 nt in length with the identical 
order of a 539-nt 5’ untranslated region (UTR), 18,788-nt open reading 
frame 1la/1b (ORF1a/1b), 3,480-nt S gene, 252-nt E gene, 654-nt M 
gene, 285-nt NS6 gene, 1,029-nt N gene, 603-nt NS7 gene and 390-nt 3’ 
UTR. The entire genome of the other three representative PDCoV strains 
(HNZK-02, HNZK-04, HNZK-06) comprised 25,419 nt in length and 
containing a 18,803-nt ORFla/1b, which had a 6-nt and 9-nt insert at 
position 1739 and 2804 nt when compared with CH-01 strain (the nt 
positions were numbered according to the entire gene of PDCoV CH-01 
sequence). The four sequenced PDCoV S genes had a 3 nucleotides 
deletion at position 52 amino acid (aa) when compared with the US 
strains and early Chinese strains (CHN-AH-2004, HKU15-44), and these 
nucleotides deletion were also existed in some Chinese PDCoV strains 
(CHJXNI2, CHN-HB-2014, CHN-HN-2014, CHN-JS-2014, HKU15-155). 
Compared with the other PDCoV strains, PDCoV CH-01 strain had 5 aa 
mutations in S protein, these changed positions were as following: His 
changed to Gln at position 99, His changed to Gln at position 148, Asp 
changed to Ala at position 260, lle changed to Leu at position 636, and 
Phe changed to Leu at position 929 (the aa positions were numbered 
according to the S gene of PDCoV CH-01 sequence). 

The PDCoV complete genome sequences examined in this study 
shared 98.4%-100% nucleotide identities with each other, and they 
shared 98.1%-99.4% nucleotide identities and 98.1% 99.4% aa 
homology with other PDCoV strains available in GenBank. The S and N 
gene of PDCoV CH-01 strain had the lowest nucleotide similarity with 
the other PDCoV strain, which shared 83.9%~95.6% and 98.4% 99.2% 
nucleotide identity with others, respectively. 


3.3. Phylogenetic analysis of PDCoV Henan strains 


To study the genetic relationships of identified PDCoV strains in this 
study with other reference PDCoV strains, phylogenetic analyses were 
first carried out based on the complete genome of the four sequenced 
PDCoV strains and other coronaviruses obtained from GenBank. The 
results showed that the sequenced Henan PDCoV strains were clustered 
into the deltacoronavirus group, forming a branch with other ten swine- 
origin PDCoV strains, which was distinct from the bird-origin delta- 
coronaviruses (Fig. 3A). Further phylogenetic analysis based on the 
complete genomes of PDCoV indicated that PDCoV CH-01 strain iso- 
lated in Henan was grouped within the Chinese clade and closely re- 
lated to CHN-HB-2014 strain (KP757891) and HKU15-44 strain 
(JQ065042). (Fig. 3B). However, the PDCoV HNZK-02, HNZK-04 and 
HNZK-06 strains were located in the same separate subcluster and more 
related to PDCoV CHJXNI2 strain (KR131621) and CH-SXD1-2015 
strain (KT021234), indicating that the origin of these isolated strains 
may be different (Fig. 3B). 

As shown in Fig. 4, phylogenetic trees were also constructed using 
the full-length PDCoV S and N genes of the sequenced Henan strains 
and other reference PDCoV strains. Phylogenetic analyses of the PDCoV 
N genes revealed that they belonged to the same groups as in the tree 
based on the PDCoV complete genomes. Phylogenetic analysis based on 
the entire S genes of PDCoV showed that CH-01 strain was closely re- 
lated with CHN-HB-2014 and HKU15-44 strains (Fig. 4A). As the 
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Fig. 1. Map of PDCoV distribution in pigs in 17 regions of Henan province. The positive rates of PDCoV infection were labeled in the regions. The different color 
scales represent different positive rates and the regions with dark colors represent high positive rates for PDCoV infection. 


phylogenetic tree based on the PDCoV complete genome, the other 3 outbreaks in China. 

sequenced PDCoV Henan strains (HNZK-02, HNZK-04 and HNZK-06) In our study, we found that the PDCoV prevalence in Henan diar- 
were still closely related each other and located in the same separate rheal pigs was high, which is similar to the other independent reports 
subcluster, and they also showed a close relationship with CHJXNI2 and from the China (Song et al., 2015; Mai et al., 2018). RT-PCR detection 
CH-SXD1-2015 strains. (Fig. 4B). Altogether, the phylogenetic data of 430 porcine diarrhea samples collected from April 2015 to March 


indicated that these sequenced strains of PDCoV in Henan were con- 2018 revealed that the total positive rate of PDCoV infection in Henan 
sistent with epidemic strains in China and might evolve from different pigs was 23.49%. Previous reports have shown that the prevalence of 
origins. PDCoV in diarrheal pigs in Guangdong and Jiangxi were 21.8% and 


33.71%, respectively (Song et al., 2015; Mai et al., 2018). In addition, 
the molecular-based detection also showed a relatively high prevalence 
of PDCovV infection in South Korea pigs and the total positive rate was 
. . ; . . _— 19.03% (Jang et al., 2017). The prevalence of PDCoV in Tibetan diar- 

Piglet diarrhea is an important disease that threatens the PIs ae rheal pigs surrounding the Qinghai-Tibet Plateau of China was sig- 
dustry development, and the coronaviruses are currently the main pa- nificantly lower than that in other provinces and the positive rate was 
thogens that caused piglet diarrhea. PDCoV is a novel swine intestine only 3.70% (Wang et al., 2018). For the limitation of the collected 
coronavirus and has caused serious economic losses for pig industry 


4. Discussion 


: : sample sizes, especially in the regions of Xu Chang (3), Luo He (3), An 
(Hu et al., 2015). So far, PDCoV has been reported in several countries Yang (6), Ji Yuan (6) and Xin Yang (6), the PDCoV existence in the pig 
and the pathogenicity for this novel virus has been studied using the herds need to be verified by evaluating more samples. So we will 
neonatal piglets (Wang et al., iia. Wane = al., 2015; Dong et al., continue to surveil the PDCoV prevalence in these regions. In addition, 
20% 6; Hu et al., 2016). However, little regarding the molecular epide- we only have collected the samples from pigs with diarrhea but not 
miology of PDCoV has been reported in Henan province, a major pig- collected the samples from pigs without such symptoms to detect 
raising region of China. In the present study, we investigated the PDCoV PDCoV infection. We will also collect the samples from pigs without 
prevalence in the diarrheal pigs in different regions and analyzed the diarrhea to fully evaluate the prevalence of PDCoV in Henan pig herds. 


full-length genome Sequence Properties of PEON. strains in Siac RT-PCRs detection indicated that the prevalence of PDCoV in the 
province, which might provide more insight about the pigs diarrhea 
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Fig. 2. Detection of PDCoV infection in pigs in Henan, China. (A) The total positive rates of PDCoV infection in diarrhea pigs from 2015 to 2018 in Henan. Two 
asterisks indicate significant differences between groups (p < 0.01 as determined by chi-square test). (B) The positive rates of PDCoV infection in diarrhea pigs at 
each season from 2015 to 2018 in Henan. Two asterisks indicate significant differences between groups (p < 0.01 as determined by chi-square test). (C) The positive 
rate of PDCoV infection in pig herds at different stages of growth. The corresponding positive rates of PDCoV infection were labeled in the histogram. (D) The positive 
rate of PDCoV co-infection with PEDV and TGEV in 101 PDCoV positive pig herds. The corresponding positive rates of co-infection were labeled in the histogram. 


diarrheal pigs was increasing year by year from 2015 to 2018, leading 
to the increase of prevalence and severity of diarrhea. Our result also 
suggested that the positive rate of PDCoV infection in spring and winter 
was striking higher than the positive rate in summer and autumn and 
the prevalence of PDCoV in the diarrheal pigs reached a high level in 
the winter, indicating that this novel virus infection mainly occurred in 
the cold season. Furthermore, the detection showed that the positive 
rate of PDCoV infection in suckling piglets was highest and the rate was 
36.43%. Interestingly, our results also showed that PDCoV can be dis- 
covered in sows and the positive rates were 2.59% in faecal samples and 
25.59% in serum samples, indicating that the PDCoV can infect sows 
and infection rate was increasing. Previous reports showed that PEDV 
can infect the pigs of all growth stages, including sows and caused se- 
vere vomiting and watery diarrhea (Yu et al., 2018). However, PDCoV 
mainly infected the piglets but rarely infected the sows in previous 
reports (Liu et al., 2018; Wang et al., 2018). Previous report suggested 
that the PDCoVs were only detected in Tibetan pigs under 1 month of 
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age but PEDV could detected in pigs of different ages (Wang et al., 
2018). A recent report found that the PDCoV can infect the sows and 
the infection rate was up to 22.4% in sows (Mai et al., 2018). In this 
study, we also proved this result, indicated that we will need to con- 
tinue monitoring the infection of PDCoV in sows. 

The co-infections of PDCoV with PEDV and TGEV in natural infec- 
tions were common in previous reports (Song et al., 2015; Jang et al., 
2017; Mai et al., 2018). Our detection also found that co-infections of 
PDCoV with PEDV were common in Henan diarrheal pigs and 61 of 101 
(60.40%) PDCoV positive samples were proved to be positive for PEDV. 
The co-infections of diarrhea coronavirus have caused severe mortality 
in piglets (Song et al., 2015) and further study will needed to elucidate 
the prevalence of these viruses in Henan province. 

Homology analysis showed that the complete genome sequence of 
PDCoV CH-01 strain shared 98.1%~99.4% nucleotide identity with all 
reference strains described in previous reports. Furthermore, the ana- 
lysis results showed that the S and N gene of PDCoV CH-01 shared 
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Fig. 3. Phylogenetic analysis based on the complete genome sequences of PDCoV strains. (A) The nucleotide sequences of the full-length genomes from four 
coronavirus genera (Alphacoronavirus, Betacoronavirus, Gammacoronavirus and Deltacoronavirus). (B). Phylogenetic analyses of complete nucleotide sequences of 4 
PDCoVs from this study and some published PDCoV sequences in Genbank. The phylogenetic tree was constructed from the aligned nucleotide sequences using the 
neighbour-joining method with MEGA 6.06 software (http://www.megasoftware.net). Bootstrap values were calculated with 1000 replicates. Reference sequences 
obtained from GenBank are indicated by strain names and GenBank accession numbers. The PDCovV isolates identified in this study are indicated with black triangles. 


Scale bars indicate nucleotide substitutions per site. 
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Fig. 4. Phylogenetic analysis using the S genes (A) and the N genes (B) from different PDCoV strains. Reference sequences obtained from GenBank are indicated by 
strain names and GenBank accession numbers. The S and N genes from PDCoV isolates identified in this study are indicated with black triangles. The phylogenetic 


tree was constructed as Fig. 3. 
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83.9%795.6% and 98.4%99.2% nucleotide identity with the reference 
strains, respectively, indicating that the S gene had a high genetic di- 
versity and strong evolutionary of PDCoV in pigs. At present, cor- 
onaviruses are divided into four genera based on the conserved do- 
mains, including Alpha-, Beta-, Gamma- and Deltacoronavirus(Chan 
et al., 2013; Ma et al., 2015). Phylogenetic analyses based on the 
complete genome showed that the Henan PDCoV strains were clustered 
into the porcine deltacoronavirus group and located in the Chinese clade. 
Further phylogenetic analysis based on the full genomes of PDCoV in- 
dicated that the CH-01 was closely related to CHN-HB-2014 strain and 
HKU15-44 strain. However, the PDCoV HNZK-02, HNZK-04 and HNZK- 
06 strains were more related to CHJXNI2 strain and CH-SXD1-2015 
strain. The results indicated that these sequenced strains might evolve 
from the different Chinese strains. The reference strains CHN-HB-2014, 
CHJXNI2 and CH-SXD1-2015 were detected from the piglets in recently 
years, indicated that these sequenced PDCoV Henan strains might 
evolve from the latest Chinese epidemic strains. In addition, the phy- 
logenetic tree based on S and N genes suggested that Henan strains 
were still clustered into Chinese clade. 

In conclusion, diarrheal samples collected from the pigs in Henan 
were detected using RT-PCR assay and evaluated for the prevalence of 
PDCoV in the diarrheal pigs. Our results indicated that the mono- 
infection of PDCoV and co-infection with PEDV in Henan pigs were 
common and the novel virus existed in pigs of different ages, including 
piglets and sows. Genetic and phylogenetic analyses suggested that the 
sequenced strains from Henan pigs were closely related to Chinese 
strains and might evolve from the latest epidemic strains in China. 
These results would provide insight into the prevalence and evolution 
of PDCoV in China and are beneficial to develop new vaccines. 
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